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In this work, we show that a quasi-one-dimensional d x 2_ y 2-w&ve superconductor with Rashba 
spin-orbit coupling is a Dili class, time-reversal invariant, topological superconductor (TS) which 
supports a Majorana Kramers Doublet (MKD) at each end of the TS. A MKD is a pair of Majorana 
end states (MESs) protected by time-reversal symmetry (TRS). An external magnetic field breaks 
TRS and drives the system from Dili to D class in which case a single MES appears at each end 
of the TS. We show that a MKD induces resonant Andreev reflection with zero bias conductance 
peak of 4e 2 /h. Experimental realizations of the proposed model are discussed. 



Introduction — A Majorana fermion is a real fermion 
which has only half the degrees of freedom of a usual 
Dirac fermion. It was first pointed out by Read and 
Green [1 that Majorana fermions exist at the vortex 
cores of 2D p x + ip y superconductors and these Majo- 
rana fermions are non-Abelian particles [2]. Soon after, 
Kitaev constructed a spinless fermion model and showed 
that a single Majorana end state (MES) exists at each 
end of a p-wave superconducting wire |3| . Recently, sev- 
eral groups proposed that effective p-wave superconduc- 
tors which support MESs can be realized when s-wave 
pairings are induced in systems with Rashba spin-orbit 
coupling [4-13^ It is predicted that these MESs induce 
resonant Andreev reflection and cause zero bias conduc- 
tance (ZBC) peaks in tunneling experiments [H|IT5]. Re- 
markably, these ZBC peaks are observed in recent experi- 
ments [T6HT8] , even though the origin of these ZBC peaks 
is still under hot debate [T9U23] . 

To realize MESs in semi-conductor wires with s-wave 
superconducting pairing, an external magnetic field is 
needed to break the Kramers degeneracy. MESs appear 
in the regime where an odd number of transverse sub- 
bands are occupied [6HT3]. In this work, we show that 
Majorana Kramers doublets (MKDs) can be realized in 
quasi-one-dimensional wires with Rashba coupling and 
d x 2_ y 2-wdive pairing, in the absence of a magnetic field. 
A MKD is a pair of MESs localized at one end of the wire 
which is protected by time-reversal symmetry (TRS) [24 . 

According to symmetry classification, without a mag- 
netic field, a quasi-one dimensional superconductor with 
d x 2_ y 2-wave pairing and Rashba terms is in Dili class as 
the system respects TRS, particle-hole symmetry (PHS) 
and breaks spin-rotation symmetry [25]. In the follow- 
ing, we show that in the topologically non-trivial regime, 
a MKD appears at each end of the wire. Interestingly, an 
external magnetic field breaks TRS and drives the sys- 
tem from Dili class to D class. In this case, a single MES 
appears at each end of the wire. The schematic pictures 
of MKDs and single MESs are depicted in Fig{l] 

It is shown previously that a single MES in a TS in- 
duces quantized ZBC peak of 2e 2 /h at zero tempera- 
ture in tunneling experiments. In this work, we show 
that a MKD in the Dili class TS induces a quantized 
ZBC peak of 4e 2 jh. We suggest that a centro-symmetric 
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FIG. 1: A quasi-one dimensional d x 2_ y 2-wave superconduc- 
tor with spin-orbit coupling, a) Majorana Kramers Doublets 
appear in the absence of an external magnetic field, b) Single 
MESs appear in the presence of a magnetic field. 



^2 -wave superconductor, CeCoIn 5 , is a candidate of 
realizing the proposed Dili class TS given that inversion 
symmetry on the surfaces is broken. 

Strictly ID model — Before studying the more real- 
istic quasi-one-dimensional quantum wires with d x 2_ y 2- 
wave pairing, we first consider a strictly one-dimensional 
version of the proposed Dili class TS. We show that the 
strictly one-dimensional model supports MKDs in the ab- 
sence of an external magnetic field. 

A Hamiltonian which describes a strictly one- 
dimensional TS and supports MKDs can be written as: 

Hid = H t + Hso + Hsc + H z 

H SO = E J>)/3 - l 2^) +1 ^{(Ty)cc^j^ + h.C. (1) 

Hsc = Ej iA (^t +1|t ^ 4 - ^UiAt) + h ' c ' 

where H tl Hso, Hsc and Hz are the kinetic, spin-orbit 
coupling, superconducting pairing and the Zeeman cou- 
pling terms of the Hamiltonian respectively. Here, ipj is 
a fermion operator at site j, a and f3 are the spin indices, 
t is the hopping amplitude, (Xr is the spin-orbit coupling 
strength, Aq is the superconducting pairing amplitude, 
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FIG. 2: a) Excitation energy versus chemical potential. The 
parameters of Hid are: L = 2000a, t — 12, Ao = 1 and 
aR = 4, where L is the length of the wire, a is the lattice 
spacing. Zero-energy modes exist in the region highlighted 
in red. b) The sum of the amplitude of the ground state 
wavef unctions, |^i| 2 + |^2| 2 , versus x where x is the site label. 
In the topologically non-trivial regime, e.g. \i — olr/2, the 
ground- state wave functions are localized at the edge, c) In 
the topologically trivial regime, e.g. \i — 1.5aR, the ground- 
state wave functions are predominantly in the bulk. 

and G y is a Pauli spin matrix. V z denotes the strength 
of the Zeeman term. It is important to note that sin- 
glet nearest neighbor pairings are introduced in Hid such 
that the pairing terms are cos k dependent in momentum 
space. This is in sharp contrast to the k independent 
s-wave (on-site superconducting) pairing introduced in 
previous works [6UT3]. 

The energy spectrum of Hid with V z = is shown 
in Figj2^i. Due to Kramers degeneracy, every state in 
Fig. 2a is doubly degenerate. It is evident from the energy 
spectrum that zero energy modes exist when the chemical 
potential satisfies \fi\ < |or|. The sum of the amplitudes 
of the two ground-state wave functions is shown in Fig|2]3 
to comfirm that the zero energy modes are end states. 
Since the ground state is doubly degenerate, we expect 
that there are two MESs, a MKD, at each end of the 
wire. In the topologically trivial regime where | /i | > |or|, 
the ground state wave- functions are predominantly in the 
bulk as shown in Fig(2]3. 

To understand the topological origin of the MKDs, we 
note that Hamiltonian Hid in momentum space can be 
written as: 

/ h(k) A(k) \ 
H ^( k >-{^(k) -h T (-k) J' (2) 

Where h(k) = {—2b cos k — /ijao + a R sm kcr y and A(fc) = 
A cos kidy . At V z = 0, the Hamiltonian respects TRS 
such that TH 1D {k)T~ 1 = H 1D {-k), and PHS such that 
PH 1D {k)P~ 1 = -H 1D {-k). Here, T = U T K and 
P = UpK, where Up = c"o ® ivy, Up = <j x (g) <jo and K 
is the complex conjugate operator. Since spin-rotation 
symmetry is also broken by the Rashba term, H\D{k) is 
in Dili class. As shown in Appendix A, H\D{k) can be 
continuously deformed into a flat band Hamiltonian Q{k) 
which is off diagonalized: 




where 

q(k) = i[e ie -( fc Vo " t v ) + e^>(<x + a y )], (4) 

and e i0±<yk ^ = -^tcos(k)-fi±a R sin(fc)+iA cos(fc) 

y/[-2tcos(k)-n±a R sin(fe)] 2 + [A cos(/c)] 2 ' 

The Dili class Hamiltonian can be classified by the Z2 
topological invariant [24j [26j [27] : 

(5) 

Here, Pf denotes the Pfaffian, T = i<j y . Ndiii can be 1 or 
— 1. The system is topologically trivial when Ndiii = 1- 
When Ndiii = — 1, the system is in the topologically 
non-trivial regime and the superconducting wire supports 
a MKD at each end of the wire as shown in Ref. [24 . For 
Hm{k), it can be shown that Ndiii = — 1 when |/i| < aR 
and Ndiii = 1 otherwise. This explains the appearance 
of the zero energy modes in Fig(2^i. It is important to 
note that finite Rashba terms are essential for the ap- 
pearance of the MKDs. These MKDs are different in 
origin from the zero energy bounded states of d x i_ y i- 
wave superconductors with no spin-orbit coupling terms 
in which case the zero energy modes are not protected 
against disorder. Moreover, without Rashba terms, the 
zero energy modes do not appear on surfaces perpendic- 
ular to the x-axis [28U30] . 

It is important to note that the two Major ana fermions 
of a MKD do not couple to each other due to TRS. 
Suppose we denote the two Major ana fermions at one 
end of the wire by 71 and 72 respectively. We have 
Tji = 72 and T72 = —71, where T is the TRS opera- 
tor. The extra minus sign in the second equation reflects 
the fact that T 2 = — 1. If the two Major ana fermions 
couple to each other, the coupling term can be written 
as ^7172. However, this coupling term breaks TRS as 
Tzcj7i72T _1 = — io;7i72. Therefore, the coupling be- 
tween the to Major ana fermions of a MKD is not allowed 
as long as TRS is preserved. This is very different from 
the case of a D-class superconducting wire in which two 
Major ana end states at one end of the wire always couple 
to each other to form a finite energy end state. 

To further verify the claim of having two MESs at each 
end of the wire, we note that Hid can be block diago- 
nalized by a unitary transformation such that 

U~ l H lD U = { H + ) , (6) 

where H± = —{2t cos k + \i ± olr sin k)a z + Ao cos ka y . It 
is interesting to note that H± respect the chiral symme- 
try <j x H±<j x = —H± such that the H± are in the AIII 
class. In the basis which diagonalize o x , H± can be off- 
diagonalized as: 

- _/ q± (k)\ 
H± -[ q l(k) J' (7) 

where q±{k) = — {2t cos k + /i) =F ap sin k + i Aq cos k. De- 
fine A±{k) = e^±W = q±{k)/\q±{k)l the AIII class 
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FIG. 3: Excitation energy as a function of chemical potential. 
The parameters of Hid are: t = 12, Ao = 1, cur — 4 and V z — 
2. Single MESs appear near the band bottom at \i « —2t, 
indicated by the blue line. Double MESs appear near the 
band center at ji w 0, indicated by the red line. The insert is 
an enlargement of the gray section which indicates that there 
are two states near zero energy at ji « 0. 



Hamiltonians can be classified by the topological invari- 
ant 



N± _ 1 r dA ± {k) 



A±{k) 



(8) 



Non-zero Naiii indicates the appearance of MESs. It 
can be shown that \N^ ITI \ = 1 when \a\ < |or| and 
WaiiiI = ® otherwise. Since both H+ and H_ are topo- 
logically non-trivial in the regime \n\ < \aji\, there are 
two MESs in the corresponding regime. This is analo- 
gous to a time-reversal invariant 2D p±ip superconductor 
which can be regarded as two copies of spinless chiral p- 
wave superconductors and each copy of the chiral p-wave 
superconductor can be classified by Chern numbers. 

Strictly ID model with finite V z — When V z is non- 
zero, TRS is broken and the Hamiltonian is no longer 
in Dili class. In this section, we show that MESs ap- 
pear even in the presence of an external magnetic field. 
The energy eigenvalues of Hid with finite V z versus the 
chemical potential are shown in Figj3J It is interesting to 
note that the zero energy modes appear in two separate 
regimes. In the regime near the band bottom, \i & —2t, 
there is a single zero energy mode in the excitation spec- 
trum which corresponds to a single MES at each end of 
the wire. 

Remarkably, two nearly zero energy modes appear near 
the middle of the band, at \i « 0, which correspond to two 
MESs at each end of the wire. It is shown in Appendix B 
that the double MESs are stable even in the presence of 
disorder. This is in sharp contrast to a one-dimensional 
TS in D class in which only a single MES at each end of 
the wire is stable [6HT3] . 

To understand the appearance of single and double 
MESs at different chemical potential, we note that the 
strictly one-dimensional system with non-zero V z is in 
the BDI class. This is because the Hamiltonian re- 
spects a time-reversal like symmetry Tbdi such that 
TBDiH 1D (k)T~^j = H 1D (-k) even though TRS is bro- 



ken. Here, Tbdi = c"o ® ctqK and T^ DI = 1. Since PHS 
is respected as before, the Hamiltonian with finite V z is 
in the BDI class. Since one-dimensional systems in BDI 
class are classified by integers, therefore, it is possible to 
have multiple stable MESs at the end of the wire [24j [25] . 

In Appdendix B, we show that Hid with finite V z 
can be classified by a topological invariant Nbdi that 
\Nbdi\ = 1 in the regime where single MESs appear, 
e.g., when /i w —2t. On the other hand, Nbdi = 2 
near fi « where two MESs appear. The condition for 
Nbdi = -1 is: 



(2t + /i) 2 <K 2 



A 2 



(9) 



This is exactly the same condition for single MESs to 
appear in the s-wave pairing case [H [5] . The conditions 
for Nbdi — 2 is: 



OCR- 



(10) 



It is interesting to note that the conditions for MKDs 
to appear, \u\ < |or|, is reproduced in Eq 10 by setting 
V z = 0. The double MESs at finite V z can be understood 
as the descendants of the MKDs. 

It is important to note that the BDI classification ap- 
plies only in the strictly ID limit, when the symmetry 
TBDiHiD(k)T^j = H\r>{—k) is respected. In quasi-one 
dimensional case, this symmetry is broken and the Hamil- 
tonian is in the D class in the presence of V z . Therefore, 
double MESs in the quasi-one dimensional case are not 
stable if TRS is broken. 

Multi-channel case — In this section, we consider the 
quasi-one dimensional limit in which multiple transverse 
sub-bands of a wire are occupied. In the quasi-one- 
dimensional case, the Hamiltonian can be written as: 

HqiD = H t + Hso + Hsc + Hz, 

Ht = E R) d, a -*W4t+d,o^R,a + h.C.) - /i^R, a ^R,c 
HSO = En,d,a,/3 "I a ^R + d,c^ * (°W? X d)^R,/3 + h.C. 

Hsc = Er ^ A o(^R + d x ,A; ~ ^R+d x ,At) - 
Hz = Er V z (^R t ^Rt - V4u^R|) • 

(11) 

Here, R denotes the lattice sites, d denotes the two 
unit vectors d x and d y which connects the nearest neigh- 
bor sites in the x and y directions respectively. This 
model is the same as the tight-binding model in Ref.[T2] 
except for the superconducting pairing terms. The pair- 
ing terms in H q \D can be written as Ao[cos(/c :E ) — cos(A^)] 
in the momentum space. Therefore, H q id describes a 
quantum wire with spin-orbit coupling and a d x 2_ y 2-w&ve 
superconducting pairing. 

The energy spectrum of H q \D with V z = and on-site 
disorder is shown in Fig(4^i. The length of the wire is 
chosen to be much larger than the superconducting co- 
herence length L ^> t/A and the width is on the order 
of the coherence length. It is evident from Fig(4^i that 
the zero energy modes appear for a wide range of chem- 
ical potential. Due to Kramers theorem, every state in 
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FIG. 4: a) Excitation energy versus chemical potential. The 
parameters in H q \r> are 14 = 0, t = 12, Ao = 1 and cxr — 
4. On-site Gaussian disorder with vairance A§ is present. 
The width of the wire is W = 6a and the length is L = 
1400a. All the states shown are doubly degenerate. Zero 
energy modes associated with MKDs appear in a wide range 
of chemical potential. The Kramers degenerate ground states 
are depicted in red. b) The parameters in b) are the same as in 
a) except V z = 2. Non-degenerate zero energy modes appear 
even when TRS is broken. The ground state is depicted in 
blue. 



Figj4^i is doubly degenerate and the zero energy modes 
are associated with MKDs at each end of the sample [24] . 
It is important to note that the zero modes are robust 
against disorder which does not break TRS. In princi- 
ple, the topological invariant of the quasi-one dimensional 
wire can be calculated using Eq.5. In the quasi-one di- 
mensional regime, q(k) in Eq.5 will be a 2N x 2N matrix 
where 2N is the number of transverse subbands of the 
wire in the normal state. 

In the presence of the V z term, time-reversal symme- 
try is broken and the H q m is in D class. The resulting 
energy spectrum of H q uj with V z = 2 is shown in Figj4j}. 
The non-degenerate zero energy modes appear for a wide 
range of chemical potential which are associated with sin- 
gle MESs at the sample end. 

Resonant Andreev reflection — It is shown previ- 
ously that a single Major ana fermion induces a quantized 
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ZBC peak of G = 2^- in Andreev reflection experiments 
[I4j [15] when a normal metal lead couples to a MES. 
Here, we show that the a MKD in Dili class TS induces 
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a ZBC peak of G = 4^- instead. 

A TS with parameters given in Fig[4] is attached to a 
semi-infinite normal lead as shown in Fig. 5a. The hop- 
ping amplitudes on the normal lead are the same as the 
hopping amplitudes on the TS. The barrier is simulated 
by a reduced hopping matrix element between the TS and 
the normal lead. Using lattice Green's function method 
[3TJ [32], the ZBC from a normal lead to a Dili class 
TS is calculated and shown in Fig. 5b. It is evident that 
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FIG. 5: a) A normal lead is attached to the end of a TS. b) 
ZBC from the normal lead to the TS versus chemical potential 
at V z =0. The ZBC is quantized at 4e 2 /h in the topologically 
non-trivial regime, due to the presence of double MESs. c) 
ZBC versus chemical potential at V z — 2Ao- The ZBC is 
quantized at 2e 2 /h in the topologically non-trivial regime due 
to the presence of a single MES. 



the ZBC is quantized at 4^- in the regime where double 
MESs appear. In the presence of an external magnetic 

2 

field, the TS is in D class and ZBC is quantized at 2^ in 
the topologically non-trivial regime as shown in Fig. 5c. 

Discussion — A few important comments follow. 
First, for simplicity, we assumed that the wire is aligned 
along the x-direction. If a quasi-one dimensional wire is 
grown along a direction tilted with angle with respect 
to the x-axis, the superconducting pairing symmetry be- 
comes cos 2#(cos k x f — cos k y t ) +sin 26 sin k x > sin k y > . Here, 
k x > and k y > denote the momenta parallel and perpendicu- 
lar to the wires respectively. In other words, one obtains 
a d x 2_ y 2 + d xy -wave pairing superconductor. Since the 
d xy term does not break the TRS and PHS, the conclu- 
sions of this work stand so long as is away from 7r/4 or 
37r/4 in which directions the pairing gap along the wire 
vanishes. 

Second, only spin singlet d x i_ y i pairing is considered 
in the main text. However, in the presence of Rashba 
terms, spin triplet pairing terms may appear [33] [34]. 
Nevertheless, spin-triplet terms do not break the TRS 
and PHS. The presence of spin-triplet terms does not 
affect the conclusion of this work as long as the bulk gap 
is not closed by these pairing terms as shown in Appendix 
C. 

Third, the results discussed in this work applies to 
all quasi-one dimensional d x 2_ y 2-wdive superconductors 
with Rashba spin-orbit coupling. A candidate mate- 
rial of Dili class TS is a layered heavy fermion super- 
conductor CeCoIn 5 . Bulk CeCoIns is a 6^.2 _^2 -wave 
superconductor [35]. Unfortunately, due to inversion sym- 
metry in the bulk, there is no Rashba spin-orbit coupling 
in the system which is crucial for the topological phases 
discussed in this paper. However, inversion symmetry is 
broken at the surface layer such that Rashba spin-orbit 
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coupling terms can be induced on the surface layer as 
shown by Maruyama et al. [36]. Therefore, the surface 
layer of a CeCoIn 5 thin film can be described by H q \D 
in Eq(TT| with the possibility of having additional triplet 
pairing terms as discussed above. It is shown in Ap- 
pendix D that multiple layers of CeCoIn 5 , with different 
Rashba strength in different layers, coupled by inter layer 
hopping terms can support MKDs. 

Finally, in Appendix F, we show that MKDs can 
emerge as end states of a metal wire with Rashba spin- 
orbit coupling if the wire is placed on top of usual 
d-wave superconductor without Rashba spin-orbit cou- 
pling. Since the d-wave superconductor is nodal, the Ma- 
jorana end state wavefunctions can leak into the d-wave 
superconductor. Fortunately, Majorana wavefunctions 
leak into the nodal directions only and the two MKDs do 
not couple to each other directly. As a result, the MKDs 
survive even if the parent superconductor is nodal. 

Conclusion — We show that quasi-one dimensional 
d x 2_ y 2-wdive superconductors with Rashba spin-orbit 
coupling are Dili class TS which support MKDs. Sin- 
gle MESs appear in the presence of a magnetic field. 
The MKDs induce resonant Andreev reflection with a 

2 

quantized ZBC peak of 4^-. We suggest that CeCoIn 5 
and metal wires on a d-wave superconductor are candi- 
date materials for this Dili topological superconducting 
phase. 
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where 



a/2 V a ° a y 



(A4) 



and D(k) = h(k) + A(k)<j y . Due to the chiral symme- 
try, the eigenvalues of the Hamiltonian can be written as 
±X a (k) with a = 1, 2. For a gapped Hamiltonian, we can 
assume A a (fc) > for all k. 

Let (xt (fe)? ? 7^(^)) T be the eigenfunctions of H\r>{k) 
with eigenvalues ±A a (fc) respectively. Using the eigen- 
value equation of Hf D (k), one obtains 

Drf X ±(k) = \lx±(k), D^Dr 1 t(k)=X 2 aV t(k). (A5) 



Therefore, the eigenfunctions of Huj(k) are 

l*a,±) = 



Xa 



where u a are the normalized eigenfunctions of DD^ . 
Once the wavefunctions are known, we can calculate the 
flat band Hamiltonian of Hir>(k), which is defined as 

Q(*0 = £a=l,2 l*a, +>(*a, +H*a, ">(*a, "|- In terms 

of u a , we have 



Q(k) 



Q(k) \_ sr ( 



t D(k) 



(A7) 



Using DDI = cos k + /i) 2 + A 2 , cos 2 k + a 2 R sin 2 k]a 
2aRsmk(2tcosk + fi)u y , we have 



q(k) 



(<t - a y ) + e ie ^ k \a + a y )}, (A8) 



Appendix A: Dili Class Topological Invariant 

In this section, we first obtain the flat band Hamilto- 
nian of Hi£>(k) in Eq.2 of the main text. We then calcu- 
late the topological invariant N^m defined in Eq.5. 

To obtain the flat band Hamiltonian, we first note that 
at V z = 0, the Hamiltonian Hir>(k) in Eq.2 of the main 
text satisfies time-reversal symmetry and particle-hole 
symmetry such that 



(-fc)c4 : 
UcH* 1D {-k)Ul 



UtH 1d 



Hm(k), and 
= -H 1D (k), 



(Al) 



where Ut = &o (g) ia y and Uc = & x (g) <jq. As a result of 
time-reversal symmetry and particle-hole symmetry, the 
Hamiltonian acquires a chiral symmetry 



UlH 1D {k)U s = -H 1D (k), 



(A2) 



where Us = iUtUc- Therefore, Hm(k) is in the Dili 
class [25]. In the basis that Us is diagonal, Hu)(k) can 
be written in the off-diagonal form 



H 1D {k) = VH 1D {k)V^ 



D(k) 
D\k) 



(A3) 



where e l6± ^ 



-2t cos(/c) — /iioijj sin(/c)+iAo cos(fc) 



= . Ac- 



y/[-2t cos(fc)-M±«R sin(ft)] 2 + [A cos(fc)] 2 ' 

cording to Refs. p^l 051 12?] . the Z2 topological invariant 
of the system can be written as: 



PfPZVA; = tt)] r 1 
^m= pf[Tg(fc = 0)] exp{-- 



dkTrtf(k)d k q(k)}}, 
(A9) 

where T = ia y . Ndiii can be 1 or —1. The sys- 
tem is topologically trivial when Ndiii = 1- When 
Ndiii = — 1, the system is in the topologically non- 
trivial regime and the superconducting wire supports a 
Majorana Kramers Doublet at each end of the wire. For 
Hio(k) in Eq.2 of the main text, it can be shown that 
Ndiii — — 1 when < and Ndiii — 1 otherwise. 



Appendix B: BDI Class Topological Invariant 

In the presence of the V z terms in Hij)(k) of the main 
text, time-reversal symmetry is broken and one cannot 
use the Z2 invariant Ndiii mentioned in the above sec- 
tion to characterize the Hamiltonian. However, we note 
that in the strictly one-dimensional case, the Hamil- 
tonian HiD(k) respects a time-reversal like symmetry 
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T B DiH 1D (k)T B ^j = H 1D (-k), where T BDI = cjq^cjqK. 
It is important to note that T^ DI = 1. Together with 
the fact that H\j)(k) respects the particle-hole symme- 
try PH 1D (k)P~ 1 = -H 1D (-k) as before, with P 2 = 1, 
H 1D (k) is in the BDI class [25]. 

It is well known that BDI class Hamiltonians in one 
dimension are classified by integer numbers [25] . In this 
section, we show how the integer topological invariant 
can be calculated following Refl371 

Due to the Tb di symmetry and the particle- hole sym- 
metry, the Hamiltonian Hm(k) acquires a chiral sym- 
metry S = a x <S> <Jo such that SHuj(k)S~ 1 = —Huj(k). 
In the basis that S is diagonal, the Hamiltonian can be 
written in the off-diagonal form 

WH 1D (k)Wt=( AT ° { _ k) A ™), (Bl) 

where 

W = 4= ( Gx ~ ax \ and (B2) 



< 
LU 




FIG. 6: Excitation energy of Hm(k) as a function of chem- 
ical potential in the presence of disorder. The parameters in 
this figure are the same as the parameters in Fig. 3 of the main 
text except the fact that on-site potential disorder with vari- 
ant Aq is added to the Hamiltonian. This is in sharp contrast 
to the multi-channel case when the system is in the D class 
in which case an even number of Majorana fermions are not 
protected against disorder. 



A(k) 



—2t cos k — ji — V z iolr sin k — Aq cos k 
-ic6R sin k + Aq cos k —2t cos k — /i + V z 



(B3) 

Note that A(k) is real at k = 0, ±7r, we can define the 
quantity 

z(k) = e i9W = Det[A(fc)]/|Det[A(jfe)]|, (B4) 

such that 0(k) = nir at k = 0, ±7r with integer n. The 
winding number of 6{k) can be used as the topological 
invariant which characterizes the Hamiltonian Huj(k). 
The winding number Nbdi can be written as 



Nbdi 



=7r dz(k) 
7T J k=0 z(k) ' 



(B5) 



It counts the number of Majorana end states at one end 
of a superconducting wire [33 [38]. Using A{k) obtained 
from Hip(k), it can be easily shown that Nbdi — 1 when 

(2t - n) 2 < V 2 - A 2 and (2t + /i) 2 > V 2 - A 2 . (B6) 
Nbdi — — 1 when 

(2t + fi) 2 < V 2 - Aq and (2t - /i) 2 > V 2 - A 2 . (B7) 

Nbdi = 2 when 

(2t ± fi) 2 > V 2 - Aq and fi 2 < V 2 + a 2 R . (B8) 

Assuming 2t^> \V Z \ and Ao, we have the Eq.9 and Eq.10 
of the main text. It is interesting to note that when 
Vz satisfies Eq |B8[ there are two Majorana end states at 
each end of the superconducting wire even when time- 
reversal symmetry is broken. These double Majorana 
end states are topologically protected and they survive 
in the presence of disorder as shown in Figj6] 



Appendix C: The effect of spin singlet and 
spin-triplet pairing terms 

In the main text, a quasi-one dimensional supercon- 
ductor with a pure d x 2_ y 2-w&ve pairing and Rashba spin- 
orbit coupling is studied. The Hamiltonian H q iD is in 
the Dili class which may support Majorana Kramers 
Doublets in the absence of an external magnetic field. 
However, due to the presence of the Rashba terms, addi- 
tional s-wave spin-singlet and p-wave single-triplet pair- 
ing channels may exist. In this section, we show that the 
proposed topological state is stable in the presence of the 
s-wave and p-wave pairing channels. 

In this section, s-wave spin-singlet and p-wave spin- 
triplet pairing terms are added to H q m in Eq.ll of the 
main text. These pairing terms can be written as: 

H As = A s £ R V4, t Vi,| + h - c - 

-MV4t+d y ,>R,a - ^R-d y ,a^R, a )] + h.C. 

(CI) 

Here, H^ s and H Ap represent the s-wave and p-wave 
pairing terms respectively, a is the spin index and 
= ±1. In the momentum space and assuming peri- 
odic boundary conditions, the pairing terms can be writ- 
ten as: 



#A 3 



A p [(sin k y + i sin k x )^\^_^ (C2) 
— (sin k y — i sin k x )iht + h.c] 



Since both of these pairing terms preserve time-reversal 
symmetry and particle-hole symmetry, adding these 
terms to H q \D of the main text does not change the 
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FIG. 7: a) Energy spectrum of a wire with d cc 2_ y 2-wave, s- 
wave, p-wave pairing and Rashba spin-orbit coupling. A s = 
A p = 0.2Ao is assumed. The parameters used are the same 
as the parameters in Fig. 4a of the main text, b) A plot of the 
sum of the amplitudes of the ground state wavefunctions of 
the wire |^| 2 = |^i| 2 + |^2| 2 in the topologically non-trivial 
regime. Here, ^i and ^2 are the Kramers pair. It is evident 
that the wavefuntions are localized at the ends c) A plot of the 
sum of the amplitudes of the ground state wavefunctions of 
the wire in the topologically trivial regime. The wavefunctions 
are predominately in the bulk. 



symmetry class of the Hamiltonian. Therefore, we ex- 
pect that the presence of Majorana end states is not af- 
fected by adding the s-wave and p-wave pairing, as long 
as these terms do not close the energy gap. In FigjT^i, 
the energy spectrum of a quasi-one dimensional wire with 
d £C 2_ 2/ 2-wave as well as the s-wave and p-wave pairings are 
shown. The ground state wave functions in the topolog- 
ically non-trivial and trivial regimes are shown in FigjTJ) 
and Fig ^jp respectively. It is evident that the Majorana 
end states are robust in the presence of the s-wave and 
p-wave pairing terms. In Figj7| A s = A p = 0.2A is 
assumed where Aq is the d x 2 _ y 2 -wave pairing amplitude. 



Appendix D: Multi-layer systems and application to 
CeCoIns 

In the main text and in the above sections, we show 
that quasi-one dimensional d x i _ y 2 -wave superconductors 
with spin-orbit coupling terms are Dili class TSs in the 
absence of an external magnetic field. In this section, 
we show that multi-layers of quasi-one dimensional d- 
wave superconductors coupled through inter-layer tun- 
neling can be a TS. In particular, we study a system with 
spatially modulated Rashba terms in which the strength 
of the Rashba terms in the top layer and bottom layer are 
non-zero but the Rashba terms in the middle layers are 
zero. We will argue below that such a model describes a 
multi-layer 6^2 _ y 2 -wave superconductor CeCoIns. 



The Hamiltonian of a multi-layer ^2 -wave su- 
perconductor with spatially modulated Rashba terms 
and inter-layer hoppings can be written as: Ht = 
Em=i H m + H tz where H m is the Hamiltonian for each 
individual layer and H tz describes inter-layer hoppings. 
Here, m is the layer label, N is the total number of layers. 
Explicitly, H m and H tz can be written as: 

Hm = H trn + HsOrn + HsCrn , 

Htm = E R ,d,c* -*(^R+d,a,m^«,m + ft ' C 

HsOm = Er^ - 29 f Bl ^R+d l a,m 4 ' (°W* X d )^R,/3,m 

+h.c. 

HsCm = E R %M(Vi +dx ,T,mV4u, m - ^R+d x ^,m^R,t,m) 

-(^R+d y ,t,m^FU,m " ^R+d y ,+,m^R,t, J + h ' C '] 

(Dl) 

Ht z = E R ,a,(m,mO ~ l z (^R,a,m^R,a,m' + h.C.) (D2) 

Here, ^R^,™ represents a fermion annihilation opera- 
tor at position R and spin a on layer m. H trni Hsom, 
and Hscm ar e the kinetic, spin-orbit coupling and the 
superconducting pairing terms respectively. It was first 
pointed out in Ref.36 that such a Hamiltonian, with the 
possibility of including small s-wave and p-wave pairing 
terms, describes multi-layers of CeCoIn 5 with spatially 
modulated Rashba terms. 

CeCoIn 5 is a layered d x 2_ y 2-wave heavy fermion su- 
perconductor. Even though many heavy fermion super- 
conductors break inversion symmetry in the bulk and 
are non-centrosymmetric superconductors, bulk CeCoIns 
respects inversion symmetry and it is not a non- 
centrosymmetric superconductor. 

However, for multi-layers of CeCoIn 5 sandwitched be- 
tween the vacuum and a substrate, the top and bottom 
layers, which are in contact with the vacuum and with 
the substrate respectively, break the mirror symmetry 
with respect to the z-axis locally. This is illustrated in 
Figj8^i. Due to the breaking of mirror symmetry with 
respect to the z-axis on the surface layers and the strong 
spin-orbit coupling of the Ce atoms, the surface layers 
acquire Rashba terms [36]. On the other hand, mirror 
symmetry of the inner layers is not broken, the inner lay- 
ers have no Rashba type spin-orbit coupling terms. As a 
result, this system has spatially modulated Rashba spin- 
orbit coupling terms [36 J. In the case of Fig{8^i, there 
are only three layers. One may assume that the Rashba 
terms in the top and bottom layers are non-zero but the 
Rashba terms of the middle layer vanishes. 

In the following, we consider a three layer system with 
^R m = (0^,0, —an/2). aR m is spatially different for dif- 
ferent layers because the difference between the vacuum 
and the substrate breaks the global inversion symmetry. 
Importantly, we consider a quasi-one dimensional geom- 
etry such that the bulk spectrum is gapped as in the 
single layer case. The energy spectrum of a finite system 
with open boundary conditions is shown in Figj8]3. It 
is evident that zero energy Majorana modes exist. The 
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FIG. 8: a) A systematic picture of CeCoIris modeled as a tri- 
layers system on a substrate. The top and bottom layers 
break the mirror symmetry with respect to the z-axis locally 
b) The energy spectrum of a finite system with open boundary 
conditions. The parameters used are Ao = 2,£ = 9Ao,aR = 
2Ao and t z = 4.5Ao c) The ground state wavefunctions in 
the topologically trivial regime at \i — —2aR. d) The ground 
state wavefunctions in the topologically non-trivial regime at 
\i — — 1.25otR. 




k In 

X 

FIG. 9: The bulk band structure for a system with peri- 
odic boundary conditions in the x-direction. There are no 
Majorana end states due to the peridic boundary conditions. 
The spectrum is fully gapped in the topologically non-trivial 
regime. The system can undergo a quantum phase transition 
from topologically non-trivial to topologically trivial phase or 
vise versa when the bulk gap is closed by tuning the chemical 
potential. 



ground state wavefunctions in the topologically trivial 
and non-trivial regimes are plotted in Fig(8]3 and Fig(8]i 
respectively. It is evident that the Majorana end states 
exist in the topologically non-trivial regime. 

In Ref.36, the authors considered a system with global 
inversion symmetry in which layers of CeCoIns are sand- 
witched between identical YbCoIn 5 layers [35 j. For ex- 
ample, in the case of a three layer system, ctR m = 
(or,0, — or) is chosen in Refl36l such that global inver- 
sion symmetry is preserved. However, it can be shown 
that such a system is topologically trivial. 



Appendix E: The Importance of the quasi-one 
dimensional geometry 

It is well known that two-dimensional d a .2_ y 2-wave su- 
perconductors are nodal and the pairing gap vanishes 
along the nodal directions due to the fact that \k x \ can 
be equal to\k y \ at the Fermi energy. However, in a quasi- 
one dimensional wire, k y is quantized and it is possible 
that \k y \ 7^ \k x \ for all k x at the Fermi energy. In this 
case, the pairing terms do not vanish at the Fermi energy 
and the system is fully gapped. 

To show that the bulk energy spectrum is gapped, we 
study a system which has periodic boundary condition 
in the x-direction and open boundary condition in the 
^-direction. The parameters are chosen to be the same 
as the finite size system in Figj4^i of the main text. Due 
to the annular geometry, there are no end states. The 
bulk excitation energy of the system versus k x at \i = 
—t is shown in Fig(9j At fi = —t, the system is in the 
topologically non-trivial regime. It is evident from Fig(9] 
that the spectrum is fully gapped. 



Appendix F: Majorana fermions in a nodal 
superconductor 

In the main text, the possibility of realizing Majorana 
fermions in intrinsic quasi-one dimensional 4 2 -y 2_wave 
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FIG. 10: The sum of the amplitude of the two ground state 
wavefunctions, \ij)\ 2 — |V ; i| 2 + IV^J 2 , of a system which consists 
of a quasi-one dimensional superconducting wire placed on 
top of a two-dimensional d cc 2_ :y 2-wave superconductor. The 
parameters are given in the text, a) The groud state wave- 
function on the quasi-one dimensional wire when the wire is 
in the topologically no n- trivial regime. It is evident that the 
ground state wavefunction is localized at the ends of the wire, 
b) The ground state wavefunctions leak into the bulk of the 
background d-wave superconductor. The wavefunctions leak 
into the nodal directions in which \k x \ — \k v \. 



superconductors are discussed. Another possible way of 
creating 4 2 - y 2 -wave pairing on a wire with Rashaba cou- 
pling is to induce d x 2_ y 2-wme superconductivity on the 
wire through proximity effect. Inducing d-wave pairing 
on wires can be experimentally challenging. In this sec- 
tion, we only discuss how Majorana fermion end states 
can survive in the presence of a nodal background. It is 
interesting to note that in this situation, the Majorana 
fermions on the wire can couple to the nodal fermions 
in the d-wave superconductor and it is not obvious that 
Majorana fermions can survive in the presence of nodal 
fermions. In this section, we show that Majorana end 
states on the wire can still survive, even though part of 



the Majorana wavefunction can leak into the d-wave su- 
perconductor. 

To show this, we couple a quasi-one dimensional wire, 
which is described by Hamiltonian H q \D in Eq.ll of the 
main text in the absence of the superconducting pair- 
ing terms, to a nodal ci^^-wave superconductor. The 
length of the wire is 300 (in units of lattice spacing) and 
the width of is 6. The d-wave superconductor has length 
500 and width 500 and the pairing amplitude on the su- 
perconductor is Ao. The parameters are chosen such that 
the dimensions of the d-wave superconductor is much 
larger than the coherence length. It is important to note 
that the pairing amplitude on the wire is zero. However, 
due to the coupling between the metal wire and the su- 
perconductor, an effective pairing is induced on the wire 
through proximity effect. The Rashba coupling strength 
on the wire is chosen to be olr = 2 Ao and on the d-wave 
superconductor. The wire is placed on top of the center 
of the d-wave superconductor. The hopping amplitudes 
on the wire and on the d-wave superconductor are chosen 
to be t = 5Ao- Each site on the wire is coupled to the 
site underneath it on the d-wave superconductor through 
hopping t z = t/2. 

The whole system is then diagonalized numerically and 
the ground state wavefunction of the whole system is 
plotted in Fi gflQ} Fig 10 1 shows the wavefunction on the 
wire and Fig|10]o shows the wavefunction on the d-wave 
superconductor. At ji = —At + 6, the wire is topolog- 
ically non-trivial. It is evident that the ground state 
wavefunction is predominantly localized at the ends of 
the wire from FigfTo|i. From FigfTo]3, one can see that 
part of the Majorana wavefunction leaks into the bulk of 
the d-wave superconductor. It is important to note that 
the wavefunction leaks into the nodal directions in which 
\k x \ = \ky\- 

Due to the presence of the gapless nodal directions, 
the Majorana fermions are no longer fully localized at 
the ends of the wire. However, since the wavefunctions 
can leak into the nodal directions only, the Majorana 
fermions at the two ends of the wire do not couple to 
each other directly in the x-direction as the x-direction is 
fully gapped. As a result, even though the energy of the 
Majorana modes is increased because of the small overlap 
of the wavefunction, this energy increase is small and the 
Majorana nature of the end states is well preserved. 
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